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We present an exact closed form expression for the finite temperature first-order density matrix 
of a harmonically trapped ideal Fermi gas in any dimension. This constitutes a much sought after 
generalization of the recent results in the literature, where exact expressions have been limited 
to quantities derived from the diagonal first-order density matrix. We compare our exact results 
with the Thomas-Fermi approximation (TFA) and demonstrate numerically that the TFA provides 
an excellent description of the first-order density matrix in the large- N limit. As an interesting 
application, we derive a closed form expression for the finite temperature Hartree-Fock exchange 
energy of a two-dimensional parabolically confined quantum dot. We numerically test this exact 
result against the 2D TF exchange functional, and comment on the applicability of the local-density 
approximation (LDA) to the exchange energy of an inhomogeneous 2D Fermi gas. 

PACS numbers: 03.75.Fi,05.30.Fk 

I. INTRODUCTION 

The recent technical advances made by DeMarco and Jin Q in the area of trapped, ultra-cold Fermi gases, has led 
to the experimental realization of what is close to being an ideal, non-interacting, many-body system of harmonically 
confined fermions. Using current state-of-the art magneto-optical traps, it is also now possible to "tune" the dimen- 
sionality of the gas from three-dimensions (3D) to quasi-2D or quasi-lD. Such a model system is of great interest to 
physicists, as it provides an opportunity to study the role of dimensionality, and the quantum statistical properties of 
a many-body system, exactly. As a result, the last few years have seen a renewed interest in the theoretical description 
of harmonically trapped ideal Fermi gases at both zero 

UM H 11 II EM and finite temperatures [H El Tlic 
primary focus of these studies has been on examining (analytically and numerically) the expressions for the local 
thermodynamic properties of the gas, e.g., the single-particle and kinetic energy densities. These quantities are, of 
course, of great importance in the density-functional theory (DFT) of inhomogeneous Fermi systems, whereby one 
can by-pass the numerically expensive one-particle Schrodinger equations. 

However, the more fundamental quantity, from which the single-particle and kinetic energy densities are both 
derived, is the first-order density matrix pi(ri,r2) 01 • Unfortunately, to date, there are relatively few examples in 
which a closed analytical form for pi(ri,i"2) can be written. One of the earliest examples dates back more than 60 
years to Husimi [l3[ . in which the zero temperature first-order densit y m atrix of a ID harmonic oscillator was derived. 
Other examples that we are aware of are the so-called Bardeen model ,14] , corresponding to a planar metal surface, and 
the work of Bhaduri and Sprung [l5| dealing with a 3D oscillator with a smeared occupancy. More recently, Howard 
et. al |16| have evaluated the zero temperature first-order density matrix of the d-dimensional harmonic oscillator for 
an arbitrary number of closed shells. Unfortunately, their form for pi(ri,r2) is somewhat impractical in that it is 
given in terms of multidimensional integrals. While these integrals are numerically easy to evaluate, they are not very 
useful for further analytical analysis (e.g, examining the asymptotic behaviour of p(r) as N — ► oo [stUHl. 

One of the central theoretical reasons for pursuing a closed form expression for /^(ri,^) is that its off-diagonal 
elements determine the exchange integrals of two-body operators, and hence, the nonlocal properties of the system. 
From an experimental point of view, a closed form expression for pi(ri,r2) is desirable because, e.g., the momentum 
density n(jp) (which is just the Fourier transform of the first-order density matrix), is experimentally accessible by 
measuring the line-shape in Compton scattering. Thus, in the case of a weakly-interacting harmonically confined Fermi 
gas, an exact knowledge of Pi(ti,T2) can serve as a benchmark from which the effects of interparticle interactions 
may be extracted. 

The primary goal of this paper then, is to present an exact closed form analytical expression for the finite temperature 
first-order density matrix of a harmonically confined ideal Fermi in any dimension. We organize our paper as follows. 
In Sec. II, we briefly review some of the basic definitions given in pUll). and then proceed to derive a closed form 
expression for the finite temperature first-order density matrix in arbitrary dimensions. Following this calculation, we 
compare our exact px(ri,T2) with the TFA in 2D and discuss the applicability of the LDA for describing the nonlocal 
properties of the 2D trapped gas. Then, in Sec. IV, we apply our results to construct a closed form expression 
for the finite temperature Hartree-Fock (HF) exchange energy density of a parabolically confined 2D quantum dot. 
We numerically investigate this exact exchange energy density and comment on its applicability in the context of the 
LDA. Finally, in Sec. V, we present a summary of our results and briefly highlight some interesting avenues for further 
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investigation. 



II. FIRST-ORDER DENSITY MATRIX IN d-DIMENSIONS 

In keeping with our our earlier work 0, 1 1 1| , we begin by considering a system of noninteracting fermions at zero 
temperature described by the time-independent Schrodinger equation 

H^{v) = [f + V(r)]Mr) = £i&(r) , (1) 

where V^(r) is a one-body potential to be specified later (all e,'s are taken to be positive). The (spinless) first- 
order density matrix can be obtained by an inverse Laplace transform of the zero-temperature Bloch density matrix, 
Co(ri,r 2 ): 



Pl ( ri ,r 2 ) = 2 J2 #(r 2 )^(ri)G(£ F ) = 

ei<E F 

where 
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-C (r 1 ,r 2 ;P) 



(2) 



C (r 1) r 2 ;/3)=^^(r 2 )^(r 1 )exp(-^ i ) , (3) 



all 



and Ep is the Fermi energy; the factor of two accounts for spin. We have put in the unit step-function <d(Ep) in 
Eq. J2J so that the Laplace transform with respect to Ep may formally be taken to be two-sided 17). Note that in 
quantum statistical mechanics, (3 is usually identified with the inverse temperature, (3 = l/fc^T. However, in our 
present context, (3 is to be interpreted as mathematical variable which in general is taken to be complex, and not the 
inverse temperature l/fc^T. 

At finite temperature, the first-order density matrix is obtained from the Bloch density matrix by using the rela- 
tion 



p 1 (r 1 ,r 2 ;T)=£; 

where 



|c T (ri,r 2 ;/3) 

P 



(4) 



CHrx.ra;/?) = C (r 1; r 2 ; /?) . ^ , (5) 

is the finite temperature Bloch density matrix, and \i is the chemical potential. In Eq. (4) , the Laplace transform with 
respect to fi is two-sided, so that fi is allowed to go negative. Specializing now to the case of an isotropic harmonic 
oscillator in d dimensions, viz., 

1 



V(r) = -mw r , r=yxj+xl-\ ^ x d > ( 6 ) 

we have for the zero-temperature Bloch density matrix |Tflj 

C'W(r 1 ,r 2 ; y 8)=cW(r+|,r-|; ( s) = (i-)^ —^jj— exp {- [r 2 tanh(/3/2) + ( S 2 /4) coth(/3/2)] } . (7) 



In the above expression (and what follows), all lengths and energies have been scaled by l osc = W%/rnu) and Tluj, 
respectively, and we have introduced the center-of-mass and relative coordinates: 

r = -(ri+r 2 ), s = r x - r 2 . (8) 

The finite temperature density matrix can, in principle, be obtained by performing the inverse Laplace transform 
given by Eq. (0J with Eq. (7). However, rather than following this direct approach (which is a very difficult task), we 
first consider the following identities 



exp{-a;tanh(/3/2)} = ^{-l) n L n {2x)e- x {e~ nP + e~ {n+1)fi } 

n=0 
oo 

exp{-ycoth(/3/2)} = ^ L k (2y)e~ v {e~ kp - e~^ +1 ^} . (9) 



k=0 
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Identifying x — r 2 and y = s 2 /4, and using (9) in Eq. (7), the Bloch density matrix now reads 

(, \ d/2 , oo oo 



OO oo 



(n+fe)/3 _ e -{n+k+2)(3 



} 



(10) 



Substituting Eq. I|1U|) into Eq. (JSJ and performing the inverse Laplace transform, Eq. |0J, leads to the finite temperature 
first-order density matrix p\ (v\ , r 2 ; T) in any dimension. For the sake of clarity, we will now proceed to give an explicit 
calculation for the simplest case of 2D, followed by a statement of the general result in arbitrary dimensions. 



A. Two dimensions 



We begin by noting the following important exact inverse Laplace transforms (all two-sided): 



C' 1 



-nf3 



sinh(/3) 



CT, 1 



2^5{j 1 -{2k + l)-n) 6(ry), 



k=0 



ttT 



sin(7r/3T) 



[exp(-f) + l] 



Putting d — 2 in Eq. I|1U|) and using Eqs. JSJ, the finite temperature first-order density matrix is given by 



P?\x,y;T) = ( ± ) £ J2(-l) n L n (2x)L k (2y)e-^ x C; 1 



-{n+k)p _ e -(n+fc+2)/3 

sinh(/3) 



TTT 



sin(7r/3T) 



(11) 
(12) 

(13) 



fc=0 n=0 

Applying the convolution theorem for Laplace transforms and making use of Eqs. (|llfl . (|12|) . we immediately obtain 



(2 
Pi 



~ ^2^2{-^) n L n {2x)L k {2y)e- {x+v) X £ / 



dr5(r - (2/ + 1) - (ra + fe)) 



1 



<1tS(t - (21 + 3) - {n + k)) 



OO OO 



^2^2(-l) n L n (2x)L k (2y)e-^ x £ 



[exp(^) + 1] 

(x+y) 



1 



[exp(^) + l] 



1 



k=0 n=0 

oo oo 



o [ [exp( 



exp( 



2l+n+k+3-fi ' 



-(x+y) 



k=Q n=0 

OG OO 



exp(^+^i) + i 



E E F^)(-irL n (2x)L k (2y)e-^ , 



(14) 



( 2) 

where the function F Tl ^ (/i) is defined as 



4%) 



exp(^^i) + 1 



(15) 



and e n — (n + 1) is the noninteracting energy spectrum (in scaled units) of the 2D harmonic oscillator potential. 
Putting back our original variables, we finally arrive at the following simple expression 



Pi 



(2)/ 



OO OO 



( r ^;T) = EE f SM^n 



x exp 



fc=0 n=0 



1- i i--2 1 2 \ r f\ri -r 2 |2 



2 |ri-r 2 | 2 



(16) 
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In terms of the center-of-mass and relative coordinates the density matrix is given by 

^( r + !' r -f ;T ) = (!)EE/S(M)(-l)'% l M^(y)exp(-(^ + ^)) (17) 

Therefore, the density-matrix only depends on the modulus of r and s, and there is a clean separation of the variables. 
Notice that Eq. (|16|l has only two sums owing to the cancellation of all but the I — term in the Z-sum of Eq. (14). In 
addition, it is readily seen that by setting r = ri = r2, we immediately obtain the finite temperature single-particle 
density 

^(n,n) ee p^(r;T) = (1) £Fi 2 )( / ,)(-l)"L„(2r 2 ) e -'- 2 , (18) 



n=0 



with 

oo 



k=Q 



[exp( 



z.-ii ' '^P( E " + T ^ ) + 1 



Equation (18) is of course identical to the result obtained in Ref. where only the diagonal part of the first-order 
density matrix was investigated. Furthermore, the exact 2D zero temperature density matrix can be obtained from 
Eq. 1)16(1 by taking the T — > limit, and when filling M + l shells, reduces to (with all dimensional constants recovered) 



M 

(2)/ 2mw v ^. . n /raw. l2 \ rl /w l2 \ / raw, 2 2 \ 

A \rur 2 ) = _53(_1)«£ B (_| ri + r 2 | 2 ) (_| ri _r 2 | 2 ) exp (- — (r? + r 2 ) 



n=0 



2mw^ / 2mw 2 \ , W 2 \ / mw / 2 J 

— i~ l^i- 1 ) L n I —f— r ) L *' ' " Vn 

71=0 ^ ' 



M-nl^ exp -— r< + - (20) 



where L]^_ n (a;) is an associated Laguerre polynomial j2£j [see also Eq. (34)] which has the property that L 1 M _ n (0) = 
M — n + 1. It is straightforward to show [with the aid of Eq. (35)] that Eq. (20) is indeed idempotent, viz., 

N= f f ^(n.raj^ra.ri) dr x dr 2 . (21) 



B. Arbitrary dimensions 

The derivation of the first-order density matrix in arbitrary dimensions closely parallels that of the 2D case. The 
w ex 
result: 



few extra steps required to obtain (r± , r 2 ; T) are clearly laid out in Ref. [llj , and so here, we simply state the final 



( 1 \ 4/2 oo oo / 1 | 2 ^ ^| | 2 

' fe = 71=0 



ri +r 2 | 2 \ r -r 2 | 2 \ _ 1 ,_ 2 2 



exp(--(rf + r 2 )) , (22) 



where 

Id) 

} ill 

v exp[(er ; + k — fx)/T] + 1 ^ exp[(er + k + 2m - M )/T] + 1 



and el d) = n + d/2. The expansion coefficients can be given in the compact form 

{d) _ 1 r(d/2 + m - 1) 



ym ml r(d/2-l) ' 1 ; 
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(2) 

In particular, we observe that gin = for all m. We have verified the correctness of Eq. (22) at T = by checking 
that it is an exact solution of the partial differential equation (valid in any dimension) 



1 <9 2 pi (£,??) 4m 2 w 2 



£77 d^d-q h 

where following the notation of [l^. we have defined £ = |ri + raj/2 and 77 = |ri — r2|/2 and 

d/2 m 



(25) 



Pi d) (£,?7) = 2 



^(-l)»i„(2£ 2 



7-1 



(M-n)/2 
m — 1 



Equation (22) also provides us with a direct way to calculate the finite temperature kinetic energy density in any 
dimension, viz., 



1. 



2(d) 



s s \ 
r + -,r -;T 

2' 2' / s=0 



(26) 



This is an easier route compared with our previous evaluation of £^ d '(r;T), which required an additional inverse 
Laplace transform of the Bloch density matrix at finite temperature (see Eq. (33) of Ref. ). It should be noted 
that when making contact with experimental studies on fully spin-polarized trapped Fermi atoms, it is appropriate 
to focus on singly filled levels, which implies that the factor of two appearing in Eq. (22) should be dropped. 

We also briefly mention that an analogous expression for p^'(ri, T2;T) for the harmonically trapped (spinless) Bose 
gas in any dimension can also be obtained by simply removing the factor of two in front of Eq. (22) and changing 
the +1 to —1 in the denominator of Eq. (23). While this resulting novel expression for the Bose gas is very different 
from what is found in the literature, it can be shown to be entirely equivalent to the more commonly used form 



Pi 



(d) 



om/T 



e -2j/T)]d/2 
2 



( , X p , - |ri + r21 tanh(j/2T) - |ri / 2| coth(j/2T) 



(27) 



C. Second-order density matrix 



In the present context of non-interacting fermions, the second-order density matrix, p^ (riT2, r^r^) can immediately 
be obtained from our knowledge of pf^ (?\ , v-i ) (this is in the spirit of the HF approximation) • An important 
quantity derived from pi is its diagonal element, which corresponds to the pair density, and is given by: 



^(nr 2 ,nr a ) = - |pW( ri )pW(r 2 ) - -\p^> {v u v 2 )\^ [closed shell] . (28) 
For example, in the HF approximation, the electron-electron potential energy is given by 

U HF f f /4^( r l r 2; r l r 2) j , / on x 

Ke = J J | ri -r 2 | dvidv2 ' (29) 

so that the the first term in (28) is responsible for the classical Coulomb repulsion and the second term yields the 
quantum-statistical exchange energy. 



III. COMPARISON WITH THE THOMAS-FERMI APPROXIMATION 



In this section, we compare our exact expression for the first-order density matrix with that of the TEA. Since the 
largest deviations between the exact and TF results are known to be at low temperatures (and small particle numbers, 
especially in low-dimensional systems), we will restrict our comparisons to zero temperature, and simply state here 
that the temperature dependence of the exact pi(r 1; r 2 ;T) is readily studied, should it be desirable. Again, for the 
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FIG. 1: Comparison of the exact zero temperature density matrix [Eq. (20) (solid curves)] with the TFA [Eq. (31) (dashed 
curves)] for various particle numbers. Panels (a)-(c) correspond to fixing the center-of-mass coordinate to r = 0, and panels 
(d)-(f) have r = 3. In all cases, pi is normalized to the central density in the trap, p(0). 



sake of simplicity, we will focus explicitly on the 2D case, although the extension of the analysis to other dimensions 
is straightforward. 

In the TFA, the 2D zero temperature Bloch density matrix is given by [l8| 



CW(ri,r 2 ;/3) 



1 

2tt/3 



e -/3y[(r 1 +r 2 )/2] e -|r 1 -r 2 | 2 /2/3 



(30) 



The inverse Laplace transform given by Eq. (2) (with Eq. (30) for the Bloch density matrix) is readily performed, 
and we find for the 2D harmonically confined gas 



2[E F -V 



ri + r 2 



Ji 



^2 (£j,-y(r))Ji y/2(E F - V(r)) 



■2[E F -V 
1 



ri + r 2 



ri -r 2 



1 



ri - r 2 



(31) 



7 



0.2 



0.0 




M = 9 




-0.05 1 — ' — ' — ' — 1 — ' — 1 — ' — 1 — ' — 1 — ' — 1 
2 4 6 8 10 12 

T 



FIG. 2: Comparison of the exact zero temperature density matrix [Eq. (20) (solid curves)] with the TFA [Eq. (31) (dashed 
curves)] for various particle numbers. Panels (a)-(c) correspond to fixing the relative coordinate to s = 1. In all cases, pi is 
normalized to the central density in the trap, p(0). Note that at the classical turning point, the TF expression drops abruptly 
to zero. 



where J\(x) is a cylindrical Bessel function and it is understood that the right-hand side of Eq. (31) is multiplied 
by the unit step function Q(Ep — V^r)). Note that Eq. (31) is also an exact solution of Eq. (25). We immediately 
observe that for Ep 3> V(r) (i.e., far from the classical turning point, r c i), Eq. (31) becomes effectively identical to 
the uniform gas result, viz., V(r) = 0. Thus, as N — > oo, r c ; — * oo, and the first-order density matrix is essentially 
unmodified by the trapping potential. 

Now, it is far from obvious by simply looking at the functional forms of Eqs. (20) and (31) that the TF and 
exact expressions are even qualitatively similar. Indeed, it is well-known that for local properties, such as p(r), the 
TF and exact densities have very different qualitative behaviour, especially in low-dimensions and small particle 
numbers. Inparticular, the shell oscillations observed in the exact expressions are not reproduced by the (smooth) 
TF profiles )4$. To address this issue, we show in Fig. 1 the exact zero temperature density matrix (solid curves) along 
with Eq. (31) (dashed curves) for various particle numbers. Let us first focus on the left panels [i.e., Fig. l(a)-(c)] 
which correspond to fixing the center-of-mass coordinate to r = 0. We observe that for low-particle numbers, there 
is a sizable quantitative difference between the two expressions, but by N ~ O(10 3 ), the exact and TF expressions 
agree quite well. In particular, the spatial oscillations in the exact and TFA expressions are in almost quantitative 
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agreement by Fig. 1(c). For the right panels (which have r = 3), we note the same trend, but now with the agreement 
between the two expressions noticeably improved (i.e., for the same number of particles). Similar behaviour is seen 
for other values of r, clearly illustrating that even at moderate particle numbers, the LDA is a good description of 
/?i(ri,r2), as a function of the relative coordinate s. 

In Fig. 2(a)-(c) we also compare the exact and TFA of p\, but now with the relative coordinate fixed to s = 1 
and the center-of-mass coordinate allowed to vary. Note that fixing s = yields the zero temperature single-particle 
density, which has already been investigated in Ref. Q. While the overall spatial behaviour of the two expressions 
are in agreement, the TFA clearly does not reproduce the fine spatial oscillations in the exact expression, which are 
associated with shell-filling effects. The reason for this, of course, is due to the fact that the LDA fails to take into 
account the discrete nature of the energy level structure of the trap. We will come back to this point in the Sec. IV 
C below. Nevertheless, by M — 49, the agreement between the exact and TF expression is quite good, except near 
the classical turning point, where the TF curves drop abruptly to zero. 

While it is possible to show analytically the reduction of the zero temperature particle and kinetic energy densities 
to their TF forms, e.g., lirn/v->oo p(r) exact — > /OTp(r) 0, E|, we have not yet been able to analytically establish a 
similar result for p(n,r2). The demonstration of this result is an interesting problem in its own right. 



IV. EXCHANGE ENERGY FOR A 2D QUANTUM DOT: HARTREE-FOCK APPROXIMATION 

As a simple application of our results, we now consider the evaluation of an exact closed form expression for the HF 
exchange energy suitable for the study of 2D parabolically confined quantum dots. Here, our motivation for focusing 
to strictly two dimensions is grounded in our previous findings for the local properties of the zero temperature trapped 
2D Fermi gas. Specifically, in Ref. 0, we showed analytically the surprising result that the 2D TF functional for the 
kinetic-energy density (i.e., without gradient corrections), when integrated over all space, leads to the exact quantum 
mechanical kinetic energy. Moreover, we also demonstrated numerically that if the exact single-particle density is 
inserted into the 2D TF kinetic-energy functional, even the local shell oscillations are reproduced remarkably well 
[2l|. Thus, the unique local and global properties of the trapped 2D system are reason enough for us to focus on 
two-dimensions. Nevertheless, we wish to re-emphasize that extending the following calculations to other dimensions 
requires nothing more than introducing the e?-dimensional measure, which is given by 

r d/2 poo 

JdM...) = Ym l ^M...). (32) 

For orientation, we will first consider the case for T — and then generalize the result to finite temperatures. 
After presenting our analytical expressions, we will close this section with some illustrative numerical results at zero 
temperature. 



A. Zero temperature 



Before proceeding with the zero temperature calculation, let us first consider a particular class of integrals that 
invariably arise during our manipulations of px{v\, r 2 ; T), irrespective of dimensionality and temperature, viz., 



(a,/3, 7 ) = / x a e- x Li(x)LZ(x) dx, (33) 



where the associated Laguerre polynomials are defined by 

_1)* 

fc! (m-k)l(J3 + k)\ 



m (— l) k (m + /3)! 



fc=0 

Inserting Eq. (34) into (33) and integrating term by term, we are left with a double sum, which can be resummed in 
closed form to give j2^| 

, „ s r(l + a)T(n + 7 + 1)T(3 - a + m) , , . 

I m , n (a, (3, 7 = > V ; 2 ( 3 F 2 l + a-/3,-n,l + a;l+ 7 ,l + a-/3-m;l , 35 

1 [m + 1)1 (n + 1)1 (1 + 7)1 [p — a) 

where 3.F2 [a, b, c; d, e; z] is the generalized hypergeometric function po|. We are now ready to proceed with the 
calculation of HF exchange energy for the 2D quantum dot. 
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The HF exchange energy, in the terms of the variables r and s, reads (hereby we set e = 1): 

1 f fh(r+|,r-|)' 2 



\\l K r + I' r -§) 



dsdr 

2 



dsdr 



where, from Eq. (20) 





2 4 



M M 



E E(- 1 )" +fei «( 2 ^ 2 )^(^ 2 )^ 2r ^-J S 2 /2)ii/- fc ( S 2 /2) e - 



2 /2 



n=0 fe=0 



Equation (36) can now be written as 

M M 



£ cx = -- E E(-l)" +fe / L n (2r 2 )L k (2r 2 )e- 2r2 dr f L\ { _ n {s 2 /2)^_ fe ( S 2 jiy^^ds 

* n=0k=0 J J ° 



(36) 



(37) 



(38) 



Going over to the variables x — s 2 /2, y — 2r 2 and making use of Eq. (35), viz. 

/„,,„(- 1/2, 1,1) = / x-^Ll^LKx^dx 
Jo 



1 1 „ 1 -, 

- — , —to, — ; 2, — n ; 1 

2 2 2 



(39) 



we obtain 



E C VL — 



— EB" 1 )^ / i„(2r 2 )i fe (2r 2 )e- 2 '- 2 dr [°° e^x-^L^xJL^xJd, 

77 n=0fc=0 ^° 

= -^EB-irtM - n + 1)(M - fc + l) r(M " " + 3/2) 

n=0 fc=0 

x / L n (y)L k (y)e- y dy 
Jo 

/^v^/w ^ 2 r(M-n + 3/2) „ 

- -^(M-n + l) 2 ^ (M _ n + ;/ s^ ; 

Equation (40) can be written in the more suggestive form 



M M 



T(M - n + 2) 



3-^2 



--,n-M, -;2,fc-M - -;1 
2 2 2 



n- M, -;2,n- M - -; 1 



E e y. 

whereby we identify the exchange energy density as 

V2 MM 



(r)dr 



(40) 



(41) 



£ x (r) = -2^V^V^(-l)«+ fe (M-n+l)(M-fc+l)L„(2r 2 )L fc (2r 2 )e- 21 ' 



n=0 fc=0 



T(M - n + ■ 



3^2 



~,n-M,^;2,k-M- 1;1 



(42) 



T(M - n + 2) 

Note that the first line in Eq. (42) is oc [^(r)] 2 . However, the terms in the second line of (42) (i.e., the hypergeometric 
and Gamma functions) prevent us from writing e x (r) as a simple functional of [/o(r)] 2 . 



B. Finite temperature 

The finite temperature exchange is readily calculated from Eq. (36) by using Eq. (16) for the finite temperature 
first-order density matrix. The calculation is entirely analogous to the zero temperature result, with the central 
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difference being that we now have to evaluate the integral [see Eq. (17)] 



7 m ,„(-l/2,0,0) = / x- 1 l 2 L m {x)L n {x)e- x dx 



r(n + l/2) 



3-F2 



r(n + l) 

Our final result for the 2D finite temperature exchange energy is given by 



1 1 1 

m,-,-;--n,l;l 



V2 



OO OC OC OO 



n=0 n'=0 fe=0 fe'=0 

x / L rl (2r 2 ) J L„,(2r 2 )e- 2r2 dr 



r(fc' + 1/2) 
r(fc' + 1) 



3F2 



fe '2'2'2 ' ' 



fc '2'2'2 ' ' 



The finite temperature exchange energy density is similarly given by 



V2 



OO OO OO CXD 



Mr;T) = -EEEEH^'^Mft^^'^) 6 " 



n=0 n'=0 fc=0 fc'=0 



r(fc' + i/2) 
x T(FTTT 3 2 



fe '2'2'2 ' ' 



(43) 



(44) 



(45) 



While Eq. (45) looks somewhat unwieldy, it turns out that the sums can be truncated relatively quickly, so that the 
temperature dependence can be numerically studied, should the need arise. 



C. Numerical results 



As mentioned above, the exact zero temperature kinetic energy density, r(r), for the trapped 2D Fermi gas agrees 
remarkably well with the corresponding TF expression, T = 7r[p(r)] 2 /2, when the exact p(r) is used as input (see 
Fig. 3 in Ref. |4|). It is then of interest to examine how well the 2D TF exchange energy density compares with the 
exact expression, Eq. (42), in both the small and lar ge N - limit. Specifically, we will investigate the applicability of 
the well-known 2D Dirac exchange energy functional 23], 

= J eJ x F (r) dr , (46) 

for the inhomogeneous 2D gas when the exact density, as given by Eq. (20) with ri = Y2, is used as input. To 
this end, we present in Fig. 3 (a)-(c) the exact (solid curves) and TF (dashed curves) exchange energy densities for 
various particle numbers. It is clear from this figure that while both expressions have similar qualitative behaviour, 
the fine spatial oscillations (i.e., shell- filling effects) are not well reproduced by the TEA. This result is entirely 
expected considering our discussion in Sec. Ill in connection to Fig. 2(a)-(c). Specifically, in obtaining e^f (r), we 
have integrated over s, so that we are left only with the center-of-mass coordinate r. Therefore, even though the 
TEA for p\ is quite good for the s integration, its failings become apparent upon an examination of the exchange 
density as a function of r; the TFA for £ x(f ) leads to a smooth, monotonically increasing profile because the resulting 
r-dependence is oc — [p TF (r)] 3//2 = —[(Ep — F(r))/7r] 3 / 2 [see Eq. (46) above]. However, in Fig. 3, ejf (r) has been 
obtained by using the exact single-particle density as input, which does have the shell effects encoded in its spatial 
dependence. Thus, although shell oscillations can be included in the TFA of e^f (r) (i.e., by using the exact density), 
it is important to note that for a given r, the weighting given to the Laguerre polynomials in [p oxact ( r )] 3 / 2 j s no t the 
same as the weight assigned to the Laguerre polynomials in Eq. (42). Consequently, the exchange energy densities 
do not agree as we vary the center-of-mass coordinate. We note in particular that the shell-effects in TFA are less 
pronounced than the exact result, and by M = 49 and r >3, are essentially washed out [see Fig. 3(c)]. It is also 
worth recounting here, that by M — 19 (see Fig. 3 in Ref. Q), the deviations between the exact and TF 2D kinetic 
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energy densities are essentially nonexistent; this is clearly not the case for the exchange energy density. The reasons 
behind the success of the TFA for the kinetic energy density are thoughoroughly discussed in Ref. Nevertheless, 
the numerics clearly indicate here that as N — > oo, £ex act ( r ) — * sJx( r )- 

While the local behaviour of the exact and TF exchange energy densities are not well reproduced (as opposed to 
the comparatively superb agreement between the exact and TF kinetic energy densities), the exchange energy itself 
agrees remarkable well. To illustrate this, we show in Table I the exact and TF exchange energies for various numbers 
of filled shells. The largest relative percentage error, which occurs at M + 1 = 10, is only AE/E ~ 0.5%. Thus, when 
global quantities are considered, the LDA is an excellent approximation for the inhomogeneous 2D exchange energy, 
even for N ~ O(10 2 ). 



V. SUMMARY AND FUTURE WORK 



The main result of this paper can be summarized by Eq. (22), which defines the exact, d-dimensional, finite 
temperature first-order density matrix of an ideal gas of harmonically confined fermions. This expression should 



12 



M + 1 


7-1 ex act 


-^CX 


AE/E 


10 


-171.71 


-170.81 


0.5% 


20 


-914.05 


-912.43 


0.2 % 


50 


-8703.06 


-8699.51 


0.04 % 



TABLE I: Comparison of the zero temperature exchange energy for the exact [Eq. (40)] and TF [Eq. (46)] expressions for 
various numbers of filled shells. The last column displays the relative percentage error in the two quantities. 

prove to be of interest in the general area of the DFT of inhomogeneous Fermi systems at both zero and finite 
temperatures. In this paper, we have used it to illustrate that (in 2D) the LDA is an excellent approximation for 
the the off-diagonal first-order density matrix in the large-A^ limit. We have also obtained a simple, closed form 
expression for the finite temperature exchange energy density for a 2D parabolically confined quantum dot. In the 
spirit of our previous work [4J , we have utilized this exact expression to test the validity of the LDA for the exchange 
energy density. In contrast to our earlier findings |4| , the 2D TF exchange energy functional does not reproduce the 
shell effects of the exact result very well. Nevertheless, when the TF exchange energy density is integrated over all 
space, we find that the resulting exchange energy is always within < 0.5% of the exact result. 

We wish to point out that the utility of our results are not limited to the topics discussed in this paper. The 
simple, analytical expressions that we have provided will be very useful in the obtaining other closed form expressions 
of interest to both theorists working in formal DFT, and experimentalists studying e.g., ultra-cold trapped fermions. 
One example that comes to mind is the 2D weakly interacting trapped Fermi gas. In Ref. the use of a contact 
psuedopotential (in 2D) led to the surprising discovery that there is no splitting between states with different angular 
momentum values I in a given shell at T = 0. This occurs in spite of the fact that the perturbation interaction does 
not preserve the SU(2) symmetry of the system. Whether this result still holds true for a finite-range psuedopotential 
is an interesting question, which can be addressed using our analytical expression for pi(ri,r2). In addition, the close 
connection between trapped Fermi gases and the theory of nuclear structure suggests that our findings will be relavent 
in the area of nuclear physics (e.g., in the Fermi gas model of the nucleus). Finally, our work here should also be 
useful in the physics of metal clusters, which represent an intermediate stage in the transition from small molecules 
to bulk solids or liquids. 
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